Optimal Condition for Strong Terahertz Radiation from Intrinsic Josephson Junctions 
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In order to enhance the radiation power in terahertz band based on the intrinsic Josephson 
junctions of Bi2Sr2CaCu2 0s+<5 single crystal, we investigate a long cylindrical sample embedded in 
a dielectric material. Tuning the dielectric constant, the radiation power has a maximum which 
is achieved when it equals the dissipation caused by Josephson plasma. This yields the optimal 
dielectric constant of wrapping material in terms of the properties of BSCCO single crystal. The 
maximal radiation power is found proportional to the product of the typical superconducting current 
squared and the typical normal resistance, or the gap energy squared divided by the typical normal 
resistance, which offers a guideline for choosing superconductor as a source of strong radiation. By 
introducing an ant i- reflect ion layer, we can build a compact device with the BSCCO cylinder and 
two wrapping dielectric layers with finite thicknesses. 

PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp 



I. INTRODUCTION 

Electromagnetic (EM) waves are vastly used in our 
daily life and their generation is considered as one of 
most developed scientific and technological fields. How- 
ever there still exists a region around terahertz (THz) 
frequencies which have wide applications such as radar, 
drug detection, safety-check, and so on, lacking of com- 
pact solid-state generator. 1,2 

It has been known for a long time that the Joseph- 
son junctions can work as oscillators to excite high- 
frequency EM waves ^ Artificial Josephson junctions were 
used first but the frequency is below THz due to 
the small superconducting energy gap in conventional 
superconductors. The discovery of intrinsic Josephson 
junctions (IJJs) in layered high-T c superconductor pro- 
vides a fantastic chance to achieve strong coherent THz 
radiationJ^ The advantages of IJJs over conventional 
low-temperature junctions are as follows. First, the junc- 
tions are homogeneous at the atomic scale guaranteed 
by the high quality of single crystals, and second the su- 
perconductivity gap is large, which in principle permits 
the frequency to cover the whole range of THz band. 
Much effort has been made to stimulate powerful radia- 
tion based on IJJs^^ - — 

An experimental breakthrough was made in 2007, 
where the THz radiation was observed based on IJJs 
from Bi2Sr 2 CaCu208+ ( 5 (BSCCO) mesa by biasing a dc 
voltage^ The key experimental results are following: the 
frequency of the radiation EM wave and the bias voltage 
obey the ac Josephson relation, and the frequency equals 
one of the cavity modes determined by the lateral size of 
mesa. 

A novel ir kink state has been proposed, which can ex- 
plain the important experimental results.— This state 
is characterized by static ±7r phase kinks in the lateral 
directions of the mesa, which align themselves alternat- 




FIG. 1: (Color online) Schematic view of proposed device with 
a thick BSCCO cylinder embedded in a dielectric material. 



ingly along the c-axis. The tt phase kinks provide a strong 
coupling between the uniform dc current and the cav- 
ity modes, which permits large supercurrent flow into 
the system at the cavity resonances, thus enhances the 
plasma oscillation and radiates strong EM wave from the 
mesa edge. 3 

Even with this novel mechanism, the radiation power 
based on IJJs is still weak due to the small radiation area 
limited by the mesa thickness (~ 1/xra), and enhance- 
ment of the radiation power remains a problem. Here we 
propose to use a long cylindrical BSCCO single crystal 
embedded in a dielectric material as shown schematically 
in Fig.l. In this way, on one hand we can increase the ra- 
diation area, and on the other hand we can have a cavity. 
Then there is a question: what is the optimal dielectric 
material for strong radiation? 
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In order to answer this question, we solve the cavity 
modes and quality factor through Maxwell's equations, 
and analyze the phase dynamics in IJJs according to the 
coupled sine- Gordon equations. We find that the radia- 
tion power has a maximum when the dielectric constant 
is tuned. It turns out that the maximum is achieved 
when the radiation power equals the plasma dissipation. 
This gives us the optimal dielectric constant of wrapping 
material in terms of the property of the superconductive 
single crystal. The maximal radiation power is propor- 
tional to the typical superconducting current squared and 
the typical normal resistance, or the gap energy squared 
divided by the typical normal resistance. 

The remaining part of the paper is organized as follows. 
In Sec. II, we discuss the cavity modes and quality factor 
of the system according to Maxwell's equations. In Sec. 
Ill, we analyze the phase dynamics associated with tt 
phase kink state in IJJs based on coupled sine-Gordon 
equations. In Sec. IV, we derive the optimal condition 
for strong radiation and the maximal radiation power 
explicitly. Sec. V is devoted to discuss the device with 
ant i- reflect ion (AR) coating. Discussion and summary 
are given in Sec. VI and Sec. VII respectively. 



II. CAVITY MODE AND Q FACTOR 
A. Cavity mode 

For simplicity we consider first a long BSCCO cylinder 
embedded in a dielectric material. The dielectric con- 
stant of wrapping material should be larger than that 
of BSCCO sample to ensure the existence of radiating 
modes. As we focus on the coherent radiation, the EM 
field is considered uniform along the c direction. 

For an open (or a radiating) system, the EM waves 
should be described by complex numbers. By consid- 
ering the Helmholtz wave equation and the continuity 
conditions of EM wave at the interface between BSCCO 
single crystal and wrapping material, the complex eigen 
wave number k rnn can be given by** 

(j^,if m ) + S r (Jm-^m) + ( 1 + &r ) j'm JmUmU'm — 0, (1) 

where Jm(y/^ck mn R) and H m (^k mn R) are the first 
kind of Bessel and Hankel functions with R the radius of 
BSCCO cylinder, e r = e/e c is the ratio between dielec- 
tric constants of wrapping material and BSCCO single 
crystal; "prime" s denote the derivatives with respect to 
R. 

There are two kinds of cavity modes, namely perfect 
magnetic conductor (PMC) like one, and perfect electric 
conductor (PEC) like one^ Here we take the PMC-like 
(1,1) mode as an example, and the extension for other 
modes is straightforward. The spatial part of electric 
field in the c direction for PMC-like (1,1) mode inside 
the BSCCO cylinder is given by^ 
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FIG. 2: (Color online) Distribution of (a) the real part and 
(b) the imaginary part of electric field in the c direction for 
PMC-like (1,1) mode. 



with the cylindrical coordinate r = (p, 0), k\\ the eigen 
wave number of PMC-like (1,1) mode, as shown in Fig. 
2. 

The angular eigen frequency, or equivalently the real 
part of eigen wave number fen in a dimensionless form, 
is given by 



UJQ 



Xiic 

RyJlT c 



(3) 



with xn ~ 1-84 standing for the first node of J[. For R = 
100/ira and e c = 16, f = uo/2tt = 0.22THz. The eigen 
frequency will not change with the dielectric constant of 
wrapping material as far as it is larger than the one of 
BSCCO sample. 



B. Quality factor 



We use quality factor 

Energy Stored 



Q = ujq 



Radiation Power 



(4) 



to characterize the energy loss caused by the radiation 
in the open system^ Considering that a cavity mode 
oscillates harmonically with time, the quality factor can 
be rewritten as 



Q = 



k^ 

2V 



(5) 



5n ( r ) = Ji(V^ fc nP) C0S< 



(2) 



where k r and k\ are the real and imaginary parts of wave 
number. The e r dependence of the quality factor can be 
evaluated from Eqs.(pQ) and (|5]). As displayed in Fig. 3, 
the quality factor increases with e r approximately in the 
form: 

Q « 0.66^-0.13. (6) 



III. PHASE DYNAMICS IN IJJS 

The phase dynamics in IJJs under bias voltage can be 
described by the inductively coupled sine- Gordon equa- 
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FIG. 3: Dependence of the quality factor on the ratio of dielec- 
tric constants between wrapping material and BSCCO single 
crystal for PMC-like (1,1) mode. 
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FIG. 4: (Color online) Dependence of current injection and 
radiation power on the bias voltage for Q=73.72 and /3 =0.05. 
The resonance takes place at V=0.45mV, which corresponds 
to the eigen frequency 0.22 THz of the PMC-like (1,1) mode 
for R=100/im. 
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tions, which are given in the dimensionless form a: 
A 7/ = (l - CA (2) ) (sin 7* + pd tli + d 2 tl i ~ Jext) , (7) 

where the lateral length is scaled by the penetration 
depth A c and time is scaled by the inverse of intrin- 
sic plasma frequency X C \^c/c; is the gauge- invariant 
phase difference at the Zth junction, ( = (\ a b/s) 2 is 
the inductive coupling which is in order of 10 5 (X a b 
the penetration depth of lateral plane and s standing 
for the period of BSCCO lattice in the c direction), 
P = 47rcr c A c / ' c^J~e~c is the normalized c-axis conductivity 
of BSCCO sample, A is the Laplace operator in lateral 
directions, A (^7/ = 7/ + i +7j_i — 27/ is the second-order 
difference operator along the c-axis. 

The general form of tt phase kink solution^!"— for the 
coupled sine- Gordon equations is given by 



Af = 2CIm(Ag) sin 7 s , 



(11) 



7 z(r, t) = ujt + A# wt + (-l)V (r), 



(8) 



where the first term is the rotating phase accounting for 
the finite dc bias voltage; the second term stands for the 
cavity term of plasma oscillation, where A is the complex 
amplitude, g(r) is the cavity mode of electric field in the 
c direction which is complex for the present open sys- 
tem as we solved above, and the frequency uo is given by 
the bias voltage based on the ac Josephson relation; the 
third term is the 7r phase kink which carries the strong 
inter j unction coupling via the I dependence. 

Inserting Eq.flSJ) into Eq.0 and omitting higher har- 
monics for I A I <C 1 we arrive at the following coupled 
equations 



Jext POJ 



/Re(Ag) cos7 s d£7 
2fdQ ' 

J g* cos7 s c^ 
2i(uj 2 -k 2 -iujf3)J\g\ 2 dn' 



(9) 



(10) 



with the integration over the lateral area of BSCCO 
cylinder. From Eq. ([TT]) it is clear that the tt phase kink 
arises at the node of cavity electric field ^ i 31 i 33 ~ — 

The radiation power and plasma dissipation per unit 
length in the c direction are given respectively by 



Pr = ^ / \E\>Xl, 



where E = Ag is the spatial part of electric field in the c 
direction. 

For a given Q factor, or equivalent ly a given dielectric 
constant (see Eq.flSJ)), the voltage dependence of current 
injection and radiation power can be evaluated by solving 
Eqs.flHI to $T2J). The result for Q = 73.73 (or e = 2.01 x 
10 5 ) is displayed in Fig. 4. 



IV. OPTIMAL RADIATION POWER 

Now we investigate the Q-dependence of the radiation 
power and the dissipation caused by plasma oscillation. 
As shown in Fig. 5, the radiation power reaches its max- 
imum when it equals the plasma dissipation. The reason 
for the existence of a maximal radiation power is that 
there exist two competing factors in the system: a good 
cavity associated with large Q makes a strong resonance 
on one hand, but it suppresses energy leak on the other 
hand. 

In order to understand the optimal condition for ra- 
diation power, we check the relations between radiation 
power and plasma dissipation. The ratio between them 
is given by 



Pdp 



(13) 
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FIG. 5: (Color online) Q-factor dependence of radiation 
power and plasma dissipation at resonance for the PMC-like 
(1,1) mode with ft — 0.05 and R — 100/ira. The radiation 
power reaches its maximum when it equals plasma dissipa- 
tion. 



where I c = J c \ 2 = c0o/87r 2 s is the typical supercon- 
ducting current with 0n the magnetic flux quantum, 
R c = l/(j c \ 2 c is the typical c-axis resistance per unit 
length of BSCCO sample. This expression looks like the 
Ohm's law but we should notice that I c is the typical 
superconducting current. 

The Eq. (|20|) can also be written as: 



Pr = 



V7T 2 A 2 

32R c s 2 e 2 ' 



(21) 



where A = ec(/>o/47r 3 A;?<7 c is the superconducting gap en- 
ergy. 

For the PMC-like (1,1) mode v « 0.15. Taking 
I c = 0.22 A, R c = 0.95 x 10 5 ^/cm, typical values for 
BSCCO, 30 the maximal radiation power is evaluated as 
85W/cm. 



From the energy conservation we have 

S 



(14) 



where S is the lateral area of BSCCO cylinder and uoj3 is 
the normal ohmic current. 

From Eqs.(|9]) and ([TO]) the super current density is given 

by 



J ext -fiu = -Re 



1LJL 



(15) 



with 



(J^cos7 5 ^) 2 + (J^cos7 5 ^) 2 

where g r and are the real and imaginary parts of cavity 
mode. From Eqs. (fT3]) to (fT5j) . the radiation power at 
resonance where uj = k r = ujq can be given by 



r 2(u /Q + f3) 2 ' 
It is then clear that P r reaches the maximal value at 

UJ 



(17) 



Q = 



(18) 



where P r P dp from Eg. ([13]) . 

By substituting Eq.(j6]) into Eq.([T8]), the optimal dielec- 
tric constant of wrapping material for strong radiation 
can be obtained 



0.20 



e c \2MnRcr c 
In Gaussian units, Eqs. (fT7|) and ([18]) give 

Pr = 77 P < • 



(19) 
(20) 



V. SETUP WITH ANTI-REFLECTION 
COATING 



In the above discussion we focus on the simple case 
in which the thickness of wrapping material is infinite, 
where the only one reflection of EM wave at the inter- 
face between BSCCO and the dielectric material deter- 
mines the cavity frequency. In what following, we show 
that two wrapping layers with finite thicknesses can work 
in the same way (Fig. 6), when the dielectric constant 
and thickness of the outside layer are tuned appropriately 
such that the two reflections at the interfaces between the 
two wrapping layers and between the outer layer and vac- 
uum cancel each other totally. The ratio between ingoing 
and outgoing of electric fields at the interface between 
two dielectric layers is given by 





, d 






i 


a 1 



BSCCO 

lelectnTmatep 
AR coating^ 



FIG. 6: (Color online) Schematic view of the proposed device 
with ant i- reflect ion coating. 
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Miii^faai)/^^ (22) 
M 2 iH£\k L a 1 )H£\k R a 2 )+M22H£ ) '(k L a 1 )H^ 

with the transfer matrix 



where /cl, kyi and /cr are the wave numbers in the inner 
wrapping layer, the AR coating and vacuum respectively, 
and d2 = ai + d (for definition see Fig. 6); "prime" s 
denote the derivatives with respect to a\ or Solving 
the equation r = the dielectric constant and thickness 
of AR coating can be determined as functions of e and a±. 
Numerical results for typical parameters are displayed in 
Fig. 7. 



VI. DISCUSSIONS 

The maximal radiation power is not affected by the 
radius of the BSCCO cylinder. The reason for this is 
that when the radius of the BSCCO cylinder increases, 
more current can be injected in, but the eigen frequency 
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FIG. 7: (Color online) Dependence of dielectric constant and 
thickness of ant i- reflect ion coating on (a) the dielectric con- 
stant of the inner wrapping layer with a± = 120/im, (b) the 
thickness of the inner wrapping layer with e — 100 of PMC- 
like (1,1) mode with eigen frequency 0.22THz. 



will decrease (proportional to the inverse of radius), with 
these two factors canceling each other. The optimal ra- 
diation power Eq. ([20]) is apparently similar to the result 
in Ref. [l]J but there the discussion was limited to non- 
resonant junctions. 

The result that the radiation power is maximal when 
it equals the dissipation can be considered as a case of 
Jacobi's law. In a recent work^i Krasnov reached numer- 
ically at a similar result for radiation in finite magnetic 
fields. This result does not depend on the geometry of 
BSCCO single crystal which can be seen from Eq.(17). 

In the present work, we focus on the case that the 
Josephson plasma is uniform along a long crystal as es- 
tablished by the cavity formed by the dielectric material. 
The maximal radiation power is therefore given by power 
per unit length. Bulaevskii and Koshelev discussed the 
case that synchronization is achieved by the coupling be- 
tween Josephson plasma and EM wave^ Therefore, their 
maximal radiation power is achieved when the length 
of BSCCO single crystal at c-axis equals the EM wave 
length. 



VII. SUMMARY 

To summarize, in order to enhance the radiation power 
in THz band we embed a long BSCCO cylinder in a di- 
electric material to increase the radiation area and mean- 
while form a cavity. We find that the radiation power has 
a maximum when the dielectric constant is tuned. The 
maximum is achieved when the radiation power equals 
the plasma dissipation which yields the optimal dielectric 
constant of wrapping material in terms of the properties 
of the superconductive single crystal. The maximal radi- 
ation is proportional to the typical superconducting cur- 
rent squared and the typical normal resistance, or the gap 
energy squared divided by the typical normal resistance, 
which offers a guideline for choosing superconductor as 
a source of strong radiation. Adding an AR coating, we 
can build a compact device with the BSCCO cylinder and 
two wrapping dielectric layers with finite thicknesses. 
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Appendix: ANTI-REFLECTION COATING FOR 
CYLINDRICAL WAVE 



For cylindrical propagation wave the electric field in 
the c direction can be written by 



E c (k,p) = A.H^ikp) - A 2 H£\kp), 



(A.1) 



where H$ and H$ are the first and the second kind 
of Hankel function, presenting the outgoing and ingoing 

I 



B e (k, p) = - UHW'ikp) - A 2 HW'(kpj\ (A.2) 



with "prime" denoting the derivative with respect to p. 
Considering the continuity of E c and Bq at the interfaces 
between the inner wrapping layer and the AR coating and 
between the AR coating and vacuum (Fig. 6), we obtain 



M^ f (k Lai ) 



A L1 
A L2 



H^'(k Mai ) 



(fc L ai) \ / 
'(fc L ai); V 

' H£\k M a 2 ) -H, 
K H^'(k M a 2 ) -H^'(k M a 2 )) 

r 



H£\k M ai) \ [A 
H^ ), {k M a 1 ) i 



AMI 
1m2 



{ m\k M a2) \ (A M i\ = ( H£\k R a 2 )\ 
L 2) 'ft M a2) VW \H^'(k R a 2 ) 



A 



Rl, 



(A.3) 
(A.4) 



where &m, and k R are the wave vectors, A^i, A\, 2 , the vacuum respectively. Eliminating Ami and Am 2 in 
and Ami, Ami, and A R1 are the amplitudes of electric Eqs.(A.3) and (A.4) we have 
field in the inner wrapping layer, the AR coating and 



M 



H^\k Lai ) 
H^'{k Lai ) 



#i 2) (fc L ai)\ (Ali 
H^'(k Lai )J V A L2 



Then we obtain the reflectance r = A^/A^i as given in 
Eq. ([22|) and the transmittance t = A^i/A^i as 



4i/7ra2 



{ n\kna 2 )\ 
M^'(k R a 2 )J 



A m . 



(A.5) 



M 21 ff<? ) (/c L ai)^ ) (fc R a 2 ) + M 22 H^ ) '(k h a 1 )H^ ) {k B ,a 2 ) - MuH™ '(k^ ai )H^ ] '' '(fc R a 2 ) - M 12 H^' '(fc L ai)ff^ ; '< \k R a 2 ) ' 



(A.6) 



where the Wronskian formula 



H&\x)H$\x)-H$Xx)Hg\x) 



-Ai 



, (A.7) 



has been used in order to obtain Eq.(A.7). 
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